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Abstract 



This dissertation reviews the Standard Model formalism as well as the Lepton 
Flavour Violating (LFV) decay processes which cause its extension, known 
as the physics beyond the SM. Firstly, using the experimental bounds on 
three body LFV decays, the corrosponding bounds on two body LFV decays 
are reviewed. The dynamical suppression of three body LFV decays due 
to momentum dependent couplings is also reviewed. Secondly, the role of 
the LFV decays to explain the LSND excess is discussed in detail, for which 
the experimental bounds on three body LFV decays, i.e. ji — > 3e are used 
to constraint the coupling gz^ e1 which is needed to calculate the anamolous 
muon decay /i — > euiUi . Then comparing the effective coupling of anamolous 
muon decay to r > 1.6 x 10~ 3 [1], it is proved that LFV is not the correct 
hypothesis to explain the LSND excess. Finally, LFV decays at loop order 
are studied in Seesaw model of neutrino masses [19] where the smallness of 
the Seesaw neutrino mass may be naturally realized with m N (mass of right- 
handed singlet neutrinos) of order 1 TeV. It is shown that the Higgs mass of 
a new scalar doublet with lepton number L — — 1 needed in the model has to 
be larger than 50 TeV to get the branching ratio of \x — > 3e to be consistant 
with the existing bound on /i — > 3e. This defeats the origional motivation of 
the model, namely that there is no physics beyond the TeV energy scale. 
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Chapter 1 
Introduction 



The most fundamental element of physics is the reduction principle. The 
large variety of macroscopic forms of matter can be traced back according 
to this principle, to a few microscopic constituents which interact by a small 
number of forces. The reduction principle has provided a guide to unraveling 
of the structure of physics from the macroscopic world through atomic and 
nuclear physics to particle physics. The laws of nature are summarized in 
the Standard Model of the particle physics. 

Standard Model is one of the successful model of the 20th century pro- 
posed by the Glashow, Weinberg and Salam. The weak and electromagnetic 
interactions are unified in the electroweak Standard Model. This model has 
provided the plenty of successful predictions with an impressive level of pre- 
cision. 

In quantum electrodynamics (QED) the interactions are specified by the 
gauge principle. The electroweak Standard Model is based on the gauge 
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symmetry group SU{2) x U(l). SU(2) is non-Abelian electroweak-isospin 
group, to which three W gauge fields are associated and U(l) is the Abelian 
hypercharge group. The associated B field and the neutral component of the 
W triplet field mix to form the photon field A and a new neutral electroweak 
field Z . This gauge group is spontaneously broken by the Higgs mechanism 
(see for example 2.3.1). In the Standard Model, all particles acquire their 
masses by interaction with another particle, the Higgs Boson. Using the 
knowledge of basic symmetry group, the gauge invariant Lagrangian has 
been written, giving not only the interactions of the various fields including 
fermions but also the mass relationships for the fermions, gauge bosons and 
Higgs boson (see for example 2.3.2). 

Despite lot of success of the Standard Model, no body can say that it is 
the end of physics. There are some limitations (see for example 2.4) which 
necessiate the extension of the Standard Model. 

When we advertize the Standard Model, we have said that it is a model 
whose foundation is symmetry. In this model only those reactions are allowed 
which conserve the individual as well as the total lepton-flavour. Thus the 
reactions violating the lepton flavour can not be contained in the Standard 
Model and cause its extension. This extension is commonly known as physics 
beyond the Standard Model. 

The problem of the physics beyond the Standard Model has been studied 
for a considerable length of time. In the past few years substantial progress 
has been made to understand this physics. Lepton Flavor Violating (LFV) 
interactions are among the most promising candidates to understand this 
new physics. 
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Using the experimental bounds on the three boby Lepton Flavor Violating 
decays we have found the bounds on the two body LFV decays. These 
bounds are suppressed by the form factors and was named as the dynamical 
suppression of the LFV bounnds. The detailed study of these LFV decays is 
presented in chapter 3. The bounds thus found on the two body LFV decays 
are used for further studies of LFV decays in some interesting reactions. 

Today there is evidence for a very important new property of the neu- 
trinos i.e. they have mass and, as a result mix with each other to lead 
the phenomenon of neutrino osscillation. Evidences that the neutrinos are 
massive particles come from three anomalous effects, the LSND excess, the 
atmospheric anomaly and solar neutrino deficit. But the atmospheric and the 
solar results are the most convincing one. The LSND has the small probability 
compared to the atmospheric and solar anomalies. In order to incorporate 
neutrino mass we must violate the Standard Model and introduce the new 
concept of LFV. A Seesaw model of neutrino masses [19] which involve right 
handed singlet neutrinos of mass m N of the order TeV, is also disscussed 
using the LFV physics at loop order. The main purpose of this model is that 
there is no new physics beyond 1 TeV. It is however, shown that even if one 
keeps ttljsj at 1 TeV, the present bounds on \i — > 3e requires that new Higgs 
boson, necessary in this model, to have mass larger than 50 TeV. This ne- 
cessiates a new physics beyond the TeV energy scale defeating the orgional 
motivation of the model. This is discussed in chapter 4, where it is also 
shown that the LFV anamolous muou decay ji — > euiUi cannot significantly 
contribute to the LSND DAR result, at least when we use g Zfie as constrained 
by /i — > 3e in a model independent way. 
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Chapter 2 

The Electroweak Theory 

2.1 History 

The proposal of the symmetry group for the Electroweak Theory, SU(2) x 
U(l), was made by Glashow in 1961. His motivation was to unify weak and 
electromagetic interactions into a symmetry group that contained U(l) em . 
The prediction includes the existance of four physical vector boson eigen- 
states, W ± , Z, and 7, obtained from the rotations of the weak eigenstate. 
In particular, the rotation by the weak angle 9 W which defines the Z weak 
boson was introduced already in this work. Furthermore the correct struc- 
ture of weak neutral current mediated by Z-boson was also obtained. The 
massive weak bosons W ± and Z were considered as mediators of weak in- 
teractions. This model has serious problem of giving masses to W ± and Z, 
since the gauge symmetry would predict their masses to be zero. The vector 
boson masses Mw and Mz were parameters introduced by hand and the in- 
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teraction Lagrangian was that of the IVB Theory. The mass term for vector 
bosons in the Lagrangian not only destroyed the gauge symmetry but also 
normalizability of the theory. 

Another key ingredient for the building of the Electroweak Theory is 
provided by the Goldstone Theorem which was introduced by Nambu in 
1960 and proved and studied with generality by Goldstone in 1961 and by 
Goldstone, Salam and Weinberg in 1962. This theorem states the existance 
of massless spinless particles as an implication of spontaneous symmetry 
breaking of global symmetries. 

The spontaneous symmetry breaking of local (gauge) symmetries, needed 
for the breaking of the electroweak symmetry £77(2) x U(l), was studied by 
P.Higgs, F. Englert and R. Brout, Guralnik, Hagen and Kibble in 1964. The 
procedure for this spontaneous breakdown of gauge symmetries is reffered to 
as the Higgs Mechanism. 

The electroweak theory as it is now known was formulated by Wein- 
berg and independently by Salam in 1968 who incorporated the gauge group 
£77(2) x U(l) introduced by the Glashow earlier. This theory, commonly 
called Glashow- Weinberg Salam Model or Standard Model (SM), was built 
with the help of the gauge principle and incorporated all the good phe- 
nomenological properities of the pregauged theories of the weak interactions, 
and in particular those of the IVB theory. It incorporated the idea of spon- 
taneous breaking of the gauge symmetry by introducing Higgs doublet. In 
this way the weak vector bosons acquaired their masses with out destroying 
the normalizability of the gauge theory.The SM is indeed a gauge theory of 
the electroweak interactions based on the gauge symmetry SU (2) x U(l) and 
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the intermediate vector bosons, 7, W and Z are the four associated gauge 
bosons. The gauge boson masses, M w and M z , are generated by the Higgs 
Mechanism in the Electroweak Theory and, as a consequence, it respects 
unitarity at all energies and is renormalizable. 

The important proof of renormlizability of gauge theories with and with- 
out spontaneous symmetry was provided by 't Hooft in 1971. 

The first firm indication that the Standard Model was the correct the- 
ory of Electroweak interaction was the discovery of weak neutral Current in 
1973 as predicted by the model. This also provided the first measurement 
of sin 2 9 W . By using this experimental input for 6 W and the values of electro- 
magnetic coupling and Gf, the SM provided the first estimates for My/ and 
Mz which were discovered experimently in 1983 at the predicted masses. 
Another important ingredients of the SM are: fermion family replication, 
quark mixing and CP violation. 

The success of the SM was clearly the discovery of the gauge bosons W ± 
and Z at the SpS collider at CERN in 1983. Since then there have been 
plenty of tests of the SM even at quantum level [2]. 



2.2 Choice of the group SU(2) x U(l) 

In order to give an argument for the choice of SU(2) x U(l) gauge group for 
electroweak unification, it is sufficient to consider the e~v component of the 
charged weak current that we write now in the form, 
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and we have introduced the lepton doublet notation and the Tj (i = 1, 2, 3) 
are Pauli matrices, 

/ 

4>L = 
Ti = 

The three generators Li of unitary, unimodular group SU (2) satisfy the com- 
mutation relation 

[7j, Ij\ i^ijklk 

i.e. the same commutation relation as satisfied by 7$ viz 




since for the fundamental representation of SU (2) 

h= T - 
2 

Note that in the charged currents there are just two generators I\ and 
I 2 . A third generator J 3 is needed in order to close the 377(2) algebera. This 
implies the third current that is relevent for the electroweak interactions, 

Obviously this can not be identified with the J e ™ which is —e^y^e. This 
clearly indicate that SU (2) is not sufficient for electroweak unification and it 



must be extended and the simple extension is to consider the group SU (2) x 
U(l). 

The Gell-Mann Nishijima relation, 

Y 

Q = h + 1 

where, 

Q = electric charge, J 3 = weak isospin, Y = weak hypercharge; implies 

Y = -1 for v L 

Y = -1 fore L 

Y = -2 for e R 

The corresponding relation amoung the current is 

rem t3 i ^ jY 

Thus 

# = 2 Or - j, 3 ) 

Therefore, if the following are used as inputs 

(-^)e L l^e L + (-1)^7^ 

, 1 s As 

{--)e L l»e L + 
(-l)e L 7 M e L + (-2)e R 7 M e R + 



jern 
J 3 



one can get Jj as output, 

Jl = 2 (j; m - 4) = 
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This clearly indicates that 

Y = -1 for 6L 

V UL 

Y = -2 for e R 

Due to symmetry breaking the two neutral currents and J^ m will mix 
to give two physical currents of which one must be identified with electro- 
magnetic current J^ m and the second current will be new currents. These 
currents will be coupled to physical vector bosons and Z^. 

9*tfW*v + = ej; m A, + gztfZp (2.1) 

where, 

= cos 6 W B^ + sin 6 w W 3fl , 

= cosOujWsh — sin 9 W B^; 6 W = weak angle. 

From here 

W 3fl = sin OyjA^ + 008 0^^, 



= cosOyjA^ - sin9 w Z fl . 



Thus, Eq. (2.1) becomes 

g 2 J* [sin 9 W A^ + cos W Z^\ + X -g Jj [cos W A^ - sin W Z^\ = eJ e ™A^ + gz^Z^. 
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So, 



e^ m = 92-^1 sin 9 W + -g cos 9 W J% = -e [e L 7 M e L + e R 7 M e R ] 
= sin W [-e L 7 M e L + ^7^] 

+ ^cos6» w [(-l)e L 7^e L + (-2)e R 7^e R + (-1)^7^] 
This simply gives us 



# 2 sin = e 
g cos 6* w = e 

Hence 
Similarly 

9z Jl = 92 J* cos Q w -\gJl sin 6 W 



tan 



52 



< / -.,,< L [ i - sin 2 0, r ^J + ^/7.-„/7. + sin 2 fl„.f /,•" r 



COS# w 

This will gives us the neutral current = Jj? c and the corresponding cou- 
pling gz- This is the main indication of the electroweak unification. 

2.3 The Electroweak Standard Model 

The Electroweak Standard Model is the commonly accepted theory of the 
fundamental electroweak interactions [3]. When we want to advertise the 
virtue of the Standard Model, we say that it is a model whoes foundation is 
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symmetry [4]. It is a gauge invariant Quantum Field Theory based on the 
symmetry group 577(2) x U(l), which is spontaneously broken by the Higgs 
mechanism. 

The Electroweak Standard Model consists of three components. 
1): The basic constituents of matter are leptons and quarks which are 
realized in three families of identical structure: 



Leptons 






V T 


e~ 


V~ 




Quarks 




u 


c 


t 




d 


s 


b 



We will concentrate on the leptonic sector only. 

2) : Four different forces act between leptons and quarks: The electro- 
magnetic and weak forces are unified in the Standard Model. The fields 
associated with these forces are spin 1 fields, describing the photon 7 and 
the electroweak gauge bosons W ± and Z. 

3) : The third component of the Standard Model is the Higgs mechanism. 
Before going to the deep discussion of the Standard Model, we have to 

require some theoretical basis [5]. 

2.3.1 The Theoretical Base 

The fundamental forces of the Electroweak Standard Model, the electromag- 
netic and the weak force, are mediated by gauge fields. The concept could 
consistantly be extended to massive gauge field by introducing the Higgs 
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mechanism which generates masses with out destroying the underlying gauge 
symmetries of the theory [6]. 

1): Gauge Sector 

Gauge invariant theories are invariant under gauge transformations of fermion 
fields: ip — > Utp. U is either a phase factor for Abelian transformations or 
Unitary matrix for non-Abelian transformation acting on the multiplets of 
the fermion field ip. Now if the theory guarantee the local transformation for 
which U depends on the space time point x, the usual space-time derivatives 
d a must be extended to covariant derivatives D which includes a new vector 
field V^. 

id^ -> iD^ = id,, - gV^ 

g defines the universal gauge coupling of the system. Under local gauge 
transformations the gauge field is transformed by a rotation plus a shift: 

V^UV^ + ig- 1 [d,U] U-\ 

But in contrast to this, the curl F of V^, 

F»„ = -ig- 1 [D^D V ] 

is just rotated under gauge transformation. 

The Lagrangian describing the systen of spin| fermions and vectorial 
gauge bosons for massless particles can be written in the compact form as 
follows: 

L V] = rjnDtl) - \ttF 2 
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It incorporates the following interactions: 
Fermion-gauge bosons 

-g$Vij). 

Three bosons couplings 

igTr (d v V^ - d^V v 
Four boson couplings 



-q 2 Tr 



2): Higgs mechanism 



What is called the Higgs mechanism is the extension of the spontaneous sym- 
metry breaking to creat massive vector bosons in a gauge invariant theory. 
As the SM is a gauge theory, the SU(2) x U(l) gauge invariance requires 
masses of the gauge bosons to be zero, since the presence of an explicit 
mass term for the gauge bosons in the Lagrangian violates gauge invari- 
ance. The Higgs mechanism circumvents this constraint by begining with a 
gauge invariat theory having massless gauge bosons. The W ± and Z° masses 
were generated by spontaneously breaking the local gauge symmetry from 
SU{2) x U(l) — > U(l) em , which was achieved by the introduction of a self- 
interacting complex scalar field, $, transforming as an SU(2) doublet. The 
doublet field <3> and its complex conjugate togather comprise four indepen- 
dent fields. Spontaneous symmetry breaking was implemented by giving one 
of the neutral fields a nonzero vacuum expectation value, 

(0>ee(O|0|O) = -|^O. 
13 



Of the four fields in the Lagrangian before spontaneous symmetry breaking, 
three fields become the longitudinal degrees of freedom of the vector bosons 
W^and Z°; the photon coupled to the remaining symmetry group U(l) em - 
generators, remains massless. 

One neutral scalar particle remains in the physical sector of the theory 
This is the so-called Salam- Weinberg Higgs particle, which the SU(2) x U(l) 
model predicts to exist. 

Since the same Higgs doublet is used to give masses to the bosons and 
fermions, which have Yukawa couplings with the scalar fields, the SU (2) x 
U(l) model predicts the couplings of the Higgs particles with all the known 
bosons and fermions but makes no prediction about its mass. This could be 
traced back to the fact that in Salam- Weinberg theory, the Higgs particle 
mass is function of the unknown quartic Higgs-boson coupling constant. 

2.3.2 Formulation of the Electroweak Standard Model 

The Matter Sector 

The fundamental fermions, as families with left handed isospin doublets and 
right handed isospin singlets appear in the fundemental representation of the 
the group 577(2) xU(l). It is realized that the symmetry pattern in the first, 
second and third generation of the fermions is same, 





VeR _ 




VpR _ 


V T 


VtR 


e~ 






L ' 







The symmetry structure cannot be derived in the Standard Model. It is an 
experimental fact that in weak interactions the parity is not conserved. The 
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different isospin assigned to the left handed and right handed field allows 
for maximal parity violation in the weak interactions. So the experimental 
observation is incorporated in the natural way. 

The relationship between the electric charge Q and basic quantum num- 
bers is described by Gell-Mann-Nishijima relation 

Y 

Q = h + T 

Interactions 

The interactions of the Standard Model are summarized by the three terms 
in the basic Lagrangian: 

L = L G + L F + L H (2.2) 
which are specified in the following way; 

Gauge fields 

SU(2) x U(l) is a non-Abelian group which is generated by the isospin op- 
erators ii, I 2 , I3 and the hypercharge Y. Each of these generlized charges is 
associated with a vector field: a triplet of vector fields W*' 2 ' 3 with J lj2j 3 and 
a singlet field with Y. The isotriplet W£, a = 1, 2, 3 and isosinglet B^ 
lead to the field strength tensors 

K» = W a + w; + ^X c 

= d^B v — dyB^ (2.3) 

g 2 is defined as the coupling constant for non-Abelian gauge group SU (2) 
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Using the equation (2.3) the pure gauge field Langragian can be written 
as follows 



L G = -\w; u W^ a -- A B^ (2.1) 



It is invariant under the non-Abelian SU{2) x U(l) transformation. 

Fermion fields and fermion-gauge interactions 

The left-handed fermion fields of each lepton family 



/ 



L 



3+ 
L 



with family index j are grouped into SU (2) doublets with component index 
a — ±, and the right-handed fields into singlets 

if = if,- 

Each left and right-handed multiplet is an eigenstate of the weak hypercharge 
Y such that the relation (2.3) is fulfilled. The Covariant derivative 

D» = dp - ig 2 IaW; + igi^B^ (2.5) 

induces the fermion-gauge field interaction via the minimal subsitution rule 

L F = J2 fti-fDntf + E jtrfD^. (2.6) 

3 3,° 

g 1 is the coupling constant for the Abelian U(l) gauge group. 
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Higgs field and Higgs interaction 



For spontaneous breaking of the SU(2) x U(l) symmetry leaving the electro- 
magnetic gauge group U(l) unbroken, a single complex scalar doublet field 
with hypercharge Y = 1 



(x) 



1 + (X) X 



(2.7) 



is coupled to the gauge fields through 



L H = (D^) + (D»<f>)-V(<f>) 



(2.8) 



B„ 



with the covariant derivative 

D^ = d lt -ig 2 IaW^ + ig 1 =f. 
The Higgs field self interaction 



\/(0) = V0 + 0+7 



(2.9) 



is constructed in such a way that has a non vanishing vacuum expectation 
value, i.e. 

1 ( 



<0> = 



V2 v 



with 



v = 



2/i 
7X' 



(2.10) 



The field (2.7) can be written in the following way, 



4(x) = 



0+ (x) 

+ H{x)+i X (x))/V2 ) 



(2.11) 
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where the field components + , H, x have vacuum expectation values zero. 

Exploiting the invariance of the Langragian, the components + , x can be 
gauged away; this means that they are unphysical (Higgs ghosts or would be 
Goldstone bosons). In this particular gauge, the unitarity gauge, the Higgs 
field has simple form 



The real field H (x) which describes small osscillations about the ground 
state defines the physical Higgs field. 

The Higgs field components have triplet and quartic self couplings follow- 
ing from V and couplings to the gauge fields via the kinetic term of Eq.(2.8). 

In addition, Yukawa couplings of the fermions are introduced in order 
to make the fermion massive. The Lagrangian for the Yukawa term can be 
written as follows 



The fermion mass terms follow from the v- part of 0° [3]. 

The Lagrangian L summarizes the laws of physics for the electromagnetic 
and the weak interactions between the leptons, and it predicts the form of self- 
interaction between the gauge fields. Morover, the specific form of the Higgs 
interaction generates the mass of the fundamental particles, the leptons, the 
gauge bosons and the Higgs boson itself, and it predicts the interactions of 
the Higgs particle [5] . 




Yukawa 



9i (v L <f> + l R + l R <\Tv L + T L <fPl R + Ir^Il) (2.12) 
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Masses and mass eigenstates of particles 

In the unitary gauge the mass terms are extracted by substiuting 



in the basic Higgs Lagrangian (2.8). The apparent 577(2) seems to be lost 
thereby, but only superficially so and remain present in the hidden form; the 
resulting Lagrangian preserves an apparent local U (1) gauge symmetry which 
is identified with the electromagnetic gauge symmetry: SU{2) x U (1) — > 



Gauge Bosons 

The mass matrix of the gauge boson in the basis Bj takes the form 



( a 2 
9w 



Ml = W 



9w 



V 



9w 9w9w 

'2 



(2.13) 



9w9w 9w ) 

This gives the mass of the vector boson in non diagonal form. The mass of 
the charged weak bosons is obvious 

M w ± = -g w v . 

As eigenstates related to the two masses M^± the charged W ± boson state 
may be defined as 

(2.14) 



W ± = — 

" y/2 



K =f K 
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For the neutral bosons, the mass term gives the matrix 



/ 



My = ~ 

v N ^ 



9w 



9w9w 

C2 



v 2 (2.15) 

\ 9w9w 9w J 

Since det(My ) = 0, therefore one of the eigenvalue of M Vn is zero. The 
above matrix is diagonalized by defining the fields A^, Z^. 



A^ = cos 6w B ^ + sin 6 wW^ 
Z^ = — sin OwB^ + cos 9 W W 3 

In matrix form the above equations can be written as follows; 

/ A 



(2.16) 
(2.17) 



Then we get 



Mi 



M 2 Z 



Where 



cos 9w sin 6 W 
— sin 6w cos 6w 



A^ : photon 

1 / 2 i '2 \ 2 



1 2 2 / 1 

4^ Uos 2 ^ 



, Q 9W 

tanfcV = — 
9w 



W 3 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



i.e. the electroweak mixing angle 9 W is defined by the ratio of the SU (2) and 
U(l) couplings. 

Introduce a parameter 



M 2 Z cos 2 e w 
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Now using the value of Mf , we get 

p=l 

This is the consquence of the fact that Higgs field is a doublet under SU(2) L 
[7]- 

Experimently the mixing angle turns out to be large, i.e. sin 2 ^ — 0.23. 
The fact that the experimental value for sin 2 Q w is far away from the limits 
0, or 1, indicate a large mixing effect. This supports the interpretation 
that the electromagnetic and the weak interactions are indeed manifestations 
of a unified electroweak interaction even though the underlying symmetry 
group SU{2) x U (1) is not simple. It may therefore be concluded that the 
electromagnetic and the weak interactions are truly unified in the Glashow- 
Salam- Weinberg theory of the electroweak interactions. 

The ground-state value of the Higgs field is releated to the Fermi coupling 
constant. From the low-energy relation Gf/V^ = Qw/^^w m P decay and 
combine with the mass relation M^ v ± = \glyV 2 , the value of v can be derived: 

v = [1/V2GF] 1 / 2 

~ 2AQGeV (2.22) 

The typical range for electroweak phenomena, defined by the weak masses 
M w and M z , is of the order lOOGeV. 
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Fermions 

The leptons are endowed with mass by means of Yukawa interactions with 
the Higgs ground state: 

M f = g f -j=. (2.23) 

Though the masses of chiral fermion fields can be introducd in a consis- 
tant way via the Higgs mechanism, the Standard Model does not provide 
predictions for the experimental values of the Yukawa couplings gf and, as 
consquence, of the mases. The theory of the masses is not available yet. 

The Higgs Bosons 

The real field H (x) which describes small osscillations about the ground 
state tells us the mass of the physical neutral scalar particles with mass 

M H = fxV2 = y/\v. (2.24) 

It can not predicted in the Standard Model since the quartric coupling A is 
an unknown parameter. 

We conclude this session with the folowing remarks: 

1) . A definite prediction of electroweak unification is the existance of 
weak neutral currents with the same effective couplings as charged currents. 
This current has been found experimentaly. 

2) . The existance of the vector bosons W ± , Z with definite masses which 
have also been discovered. 

3) . The theory has one free parameter sm 2 9 w . 
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2.4 Limitations of the Standard Model 



There is no confirmed experimental evidence from accelerators against the 
Standard Model, and several possible extensions have been ruled out. Never- 
theless, there is no thinking physicst could imagine that the Standard Model 
is the end of physics. Even if one accepts the strange set of group repre- 
sentations and hypercharge that it requires, the Standard Model contains at 
least 19 parameters. Moreover, many more parameters are required if one 
wishes to accommodate non-accelerator observations. For example, the neu- 
trino masses and mixing introduce at least 7 parameters: 3 masses, 3 mixing 
angles and 1 CP violating phase [8] . 



2.5 Summary 

In this chapter, the SU(2) and U(l) gauge group have been reviewed, which 
provides a back ground knowledge for the understanding of the Electroweak 
Standard Model. The concept of how vector bosons W ± , Z acquire masses 
have also been discussed in this chapter. The neutrino masses and mixings 
which are the weaknesses of the Standard Model have been mentioned here 
and will be studied in detail in the last chapter. 
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Chapter 3 

Generic Feature of Lepton 
Flavour Violation 

3.1 Introduction 

It is generally believed that the Standard Model of electroweak interactions 
is a low-energy approximation to a more fundamental theory. Yet there is 
no clear experimental evidence either to guide its extension to additional 
physical processes or to predict the model parameters. The Standard Model 
incorporates the lepton family-number conservation, which has been empir- 
ically verified to high precision but is not a consequence of a known gauge 
theory . Indeed many theoretical extensions to the Standard Model allow lep- 
ton family-number violation within a range that can be tested by experiment 

[!]■ 

The predications of the rate for a given family-number nonconserving 
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process vary among these extensions, and thus provide a test of the model. 
Many possibilities have been explored, and the present experimental limits 
for a wide variety of the processes have been tabulated. Thus the Lepton 
Flavour Violation (LFV) processes that are strongly suppressed in the Stan- 
dard Model by powers of (small) neutrino masses may provide signals for 
new physics. At present we have strigent bound for \i decays, e.g. 

BR[p -»• 3e) < 1(T 12 

and some what weaker O(10~ 16 ) bounds on LFV r decays [2]. 

The possibility of large samples of decaying vector bosons [V = J/ip, T, 
and Z°] and the clear signature provided by \x k e ± final states suggest search- 
ing for LFV two-body decays 

V -> /2 ± e ± . 

Here we show that rather simple considerations based on unitarity, provide 
rather strong constraints on two-body LFV processes. Hence, most three- 
body /j, and r LFV decays are likely to provide more sensitive tests of lepton 
flavour violation, rather than the two-body decays. 

3.2 Basic consideration and calculation 

Let us assume that a vector boson V* (here Vi could be either a fundamental 
state, such as the Z° or the quark- antiquark bound state such as the 0, J/ip, 
or T) couples to /i ± e ± . If it couples also to e + e~, as all the states above 
do, then by unitarity V exchange contributes also to /i — > 3e. Let us define 
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the effctive couplings between vector boson Vi and ji ± e ± as gv , and the 
corresponding effective Lagrangian can be written as 

£eff = ~9v^l a eV a + H.c. (3.1) 

From the Feynman diagram, the amplitude A (/i — > 3e) can be written as 

A(fi^3e) = u e (k 3 )gv lie 'fu li (p) j^rz s ^(h)9v ee l l3 v e (k 1 ) 

= u e (k^ a u^(p)u e (k 2 )^ a v e (ki) g ^f^\ (3.2) 

Here gv ee is the effective coupling of the vector boson Vi to e + e~, while 

s = (k 1 + k 2 f < mj. 

(There are, of course, also potential axial vector couplings of V to e + e~, 
which contribute to this process. These can be included in the above, but as 
we shall see they do not change qualitative discussion). 

Since we are dealing with the low energy process, as a first approximation, 
it is sensible to neglect s in comparision with My. Therefore the Eq. (3.2) 
takes the form, 

A(fJL — 3e) = Ueih^U^Ueih^Veih)^^. (3.3) 

2 

In order to calculate the decay width we have to calculate \A\ which is given 



by 



\A\ 2 = AA ] 

= u e {kz)^ a u^{p)u e {k 2 )i a v e {ki)v e {ki) 

x lp u e (k 2 )u,(p)Yu e (k 3 ) fe^V ■ (3.4) 



My 
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Taking the spin average for the initial state and summation over spins for 
final state particles, we get 

H l^l 2 = Ye[ Tr ha(M -m e )-f p { / k 2 + m e )} 



spins=±^ 

xTr 



2 

/ ~ \ 2 

9V ue 9V e 



+ mJYih + rn*) ( ^« ) ■ (3.5) 



Ml 



Assuming that the external particles (electrons) are massless. Therefore, Eq. 
(3.5) takes the form 

£ | A\ 2 = 2 [(h ■ p) (k 2 • fc 3 ) + (h ■ h) (k 2 ■ p)} . (3.6) 

spins=±^ 

We next carry the phase space integrations, starting with the integrals 
over the electrons momenta, given by 

(q) = I d'k^ki^S 4 (h + k 2 -q) (3.7) 

where 

q = p-h. (3.8) 

It follows from the Lorentz covariance of the integral (3.7) that the most 
general form is 

Jf" ( q ) = g^A (g 2 ) + q"q v B (g 2 ) . (3.9) 
From this equation it follows that 

g^{q) = AA(q 2 )+q 2 B(q 2 ) 

q,qJ^{q) = q 2 A (q 2 ) + (q 2 )" B (q 2 ) . (3.10) 
This can be obtained by using Eq. (3.9). 
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Since we have been taking the electron masses to be zero so that k\ = 
k 2 = k 2 = and, on account of the 5-function in (3.7), 

q 2 = (h + k 2 f = 2 (fci • k 2 ) 

\ = h-k 2 (3.11) 

In order to find A (q 2 ) and B (q 2 ), we calculate the expersion on the left 
hand sides of (3.10). From Eq. (3.7) we obtain 

9,J^{q) = ld^A 9 -^^5\k l + k 2 ~q) 

= fd^d^-^-S'ih + h-q). (3.12) 

Using Eq. (3.11), the above equation takes the form 

- i 9 2 /(V). (3.13) 

We see from its defination that the integral / (q 2 ) is an invariant, so that 
it can be evaluated in any coordinate system. For our convienence we shall 
choose the centre-of-mass system of two electrons. In this system 

ki = -k 2 

And for the massless electron 



E 1 = V|ki|+m e 2 
= |ki| 

also 

E 2 = |k 2 | . 
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So from above two results, we can write 



E 1 = E 2 = E. (3.14) 
Hence by removing the integration on fc 2 , we get 

r E 2 dE 
= j -^S(2E-q )dn 

= 2tt (3.15) 

and from Eq. (3.13) 

(q) = irq 2 . (3.16) 
Similarly, calculating q IJi q u I IJ ' 1/ (q), one gets 

j3i, j3i„ WvKK. r4 



J^(q) = fd^d^^pfldHh + h-q) 

= J d^k 2 {q - h) 5Hh + h - q) . (3.17) 
This is obtained by using Eq. (3.7). Also using Eq. (3.11), we got 

q,qJ^(q) = ^-1 (q 2 ) 

= \-{q 2 ) 2 - (3.18) 

Substituting Eqs. (3.16) and (3.18) in Eq. (3.7), we get 

-(g 2 ) = (g 2 ) + q 2 B (g 2 ) (3.19) 

\*(q 2 ) 2 = q 2 A(q 2 ) + (q 2 ) 2 B(q 2 ). (3.20) 
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Solving (3.19) and (3.20) for the vaules of A (q 2 ) and B (q 2 ), we have 

Hi 2 ) 



6 



2tt 
~6~' 



Hence Eq. (3.9) leads to 



6 

Then the partial muon decay rate can be written as 



(3.21) 



dT = 



1 1 d 3 k s r^kx^ka 



(2vr) y 2^ E 3 

2 



r d J ki 
J E\ 



Eo 



]T \A\ 2 (27r) 4 5 4 (k 1 + k 2 -q) 



spins=±i 



9V^9Ve, 
x My 



1 d 3 k a 



(2tt) & 2E m E 3 



x2[ PlI k 3 J^(q) + k 3lI p u I^ (q)} 
Using (3.21), (3.22) becomes 



(3.22) 



dT 



9v„ e 9v ei 



1 1 d 3 k 3 



M 2 



x2 



9v^gv e 



(2n) 5 2E ll E 3 
P»k 3u {g^q 2 + 2q"q v ) 
+h„Pv (9" u q 2 + 2q»q") 
°' 1 1 d 3 k 3 n 



M 2 



(2n) 5 2E ll E 3 6 



x4 



(3.23) 



(P • fca) g 2 + 2 (fc 3 • q) (p ■ q) ■ 

This is obtained by using Eq. (3.7). 

Finally, we must integrate Eq. (3.23) over all momenta k 3 of the emitted 
electron. For a muon at rest, i.e. in the rest frame of muon, we have 

p = (^,o) = K,o) 

k 3 = (E 3 ,k 3 )- (3.24) 
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Also 



go = nfy - E 3 

q = -k 3 

Now, we will calculate the terms involving in Eq. (3.23), i.e. 

q 2 = (p + hf 



(3.25) 





(ml- 






k\ — 





k 3 .p = 




h-q = 


m»E 3 


p.q = 


m M (m 



(3.26) 



By substuting these values Eq. (3.23) takes the form 



dY = 



'9v^gv ee \ 2 _L 1 ^ 3 k 3 7r 

/ (27r) 5 2m M £ 3 6^ 



x 



[m^Ez) (ml - 2m II E 3 ) + 2 (m^E^) m M (m M - E 3 ) 



, Ml 



1 2tt 
(27f"6~ 



EjdE 3 dVt 



3ml - 4m M £ 3 



(3.27) 



Integrating Eq. (3.27) over all directions Q of the emitted electron and over 
its complete range of energies < E 3 < we obtain the total decay rated 



r = 



47T 27T 



(2tt) 5 6 m "Jo 



'gv^gvA 2 4tt 27rm5 



£f [3m M - 4£ 3 ] d£ 3 



r (fj, -> 3e) 



\ M{r J (2tt) & 6 16 

i fgv ue gv e A 2 m l 



Ml 



1927T 3 
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(3.28) 



Comparing the above contribution to the /i — > 3e process to that of 
the ordionary muon decay, \i — >■ ez/z/, which proceeds via exchange and 
(almost) identical kinematics, gives the relation 

rl/i-^ew/] ~ My 'm^ 
This ratio is defined as the branching ratio. Therefore 

3e)] y _ e ^ (3.30) 



Since 



and 



while 



T(v^e+e-) ~g Ve M v , 



r -> ez/) ~ g 2 w M w , 
we can rewrite the above expression as 

[B (// -> 3e)] y _ exc ^ w — - . (3.31) 



T 2 (W -> ei/) V My 

Using [BR(fi — > 3e)]y_ e;£ . ch < 1CT 12 and other data pertaining to the e + e _ 
widths of the various vector mesons Vi, we find a set of bounds for the two 
body LFV branching ratios of these vector mesons, i.e. 
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These bounds are calculated by dividing the above decay width by the full 
width of the vector mesons. Using the numerical values these bounds become 



B(Z^fie) < 5x10 



-13 



(3.33) 



B (J/ip -> lie) < 4 x 10 



-13 



(3.34) 



B(T^fie) < 2x10 



-9 



(3.35) 



B($ -> lie) < 4 x 10 



-17 



(3.36) 



All the vector (or psudoscalars) used as intermediaries in deriving the 
bounds in Eqs. (3.33-3.36) are not on-shell. Thus we must entertain the 
possibilty that their contribution to the three-body decays considerd are 
reduced. This could at least weaken the various strong bounds obtained 
above. Now we will focus on the possible mechanism for the reduction of 
such bounds. 

3.2.1 Dynamical suppression of the LFV bounds 

The source of suppression is connected to possible "form factor" effect due to 
the dynamics which would, for example, reduce the contribution of various 
Vi states to /i — > 3e compared to the naive expectations. However the effect 
of these form factors should be minimal or controllable if the LFV is induced 
by physics at scales much higher then the EW scale of the Z mass. These 
effects of dynamics are beutifully discussed by Illuna, Jack and Rienmann 
in their recent papers [3, 4]. Following these two papers, first of all we will 
discuss the dynamical suppression in the two body decay of Z — > if and 
then move to ji — > 3e. 
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The rare process fj, — > is the classic example of the lepton flavour 
violating process. The previous limit on its branching ratio is [5] 

erf) < 4.9 x 10" 11 . 

This reaction has not been observed so far, and the best experimental upper 
limit of its branching fraction is [4] 



B(n -> e 7 ) < 1.2 x 10 



-li 



i.e. it is improved by the factor of 4.1. 

At the Z factory LEP, searches for similar processes, through the Z boson, 
became possible: 

Z — > e/i. 

The best experimental limit on its branching ratio is (95% C.L.) 



B(Z^eV) < 1.7 x 1(T 6 . 



Let's discuss this reaction in detail. 

The most general matrix element for the interaction of an on-shell vector 
boson with a fermionic current, as shown in fig. 3.2.1, may be described 
by four dimensionless form factors. At one loop level it is convinent to 
parameterize the amplitude as 

A = --^e^ e (p 2 )T p v,( Pl ), (3.37) 

2 

with aw = (f^y, £ being the boson polarization vector and 

q" 

T P = 7 P (fv - fAlb) + YF~ Mm + /b7s) <V, (3.38) 

M w 
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where 



In Eq. (3.38), the form factors fy and Ja stands for vector and axial- vector 
couplings and fu and for magnetic and electric dipole moments/transitions 
of equal/unlike final fermions. The form factors depend on the momentum 
transfer squared Q 2 = (p 2 + Pi) 2 ■ In principle all the four form factors must 
be non-zero. In order to find the decay width we have to find \A\ 2 which is 



\A\ 2 = AAl 



(3.39) 



So 



fgaw\ 
V 4tt J 



+ 



8(P2-Pl){\fv\ 2 +\fA\ 2 } 

{\f E \ 2 + \f M \ 2 } 



Ml 



{16 (p 2 • q) (pi • q) - 4 (q ■ q) (p 2 • p ± )} 



(3.40) 



Using the spin averages, above result becomes 



spins 



\A\ 2 = 




) 21 f 




V 4tt , 


12 . 



+ 



{|^I 2 + |/m| 2 } 



{16 (pa • g) (pi • g) - 4 (g • g) (p 2 • Pi)} 

(3.41) 



The decay width can be calculated as 

1 r r d 3 px d 3 p 2 

(2tt) 3 E 1 (2n) 3 E 2 



r = 



2M Z 

„3 



]T |A| 2 (2vr) 4 5 4 (g- Pl -p 2 ) 



spins 



a 



24vr 2 



Jvl 2 + |/^| 2 + tjV (l/^l 2 + l/^| 2 ) 



(3.42) 
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The branching ratio can be obtained by dividing the decay width obtained 
in Eq. (3.42) by the total width of Z boson. 



n \ 7 i & w M z 



fv\ 2 + \fA\ 2 + ^r(\f E \ 2 + \fM\ 2 ) 



(3.43) 



where 



c w 



aw 
Mw 



M z 



These form factors are model dependent. It is easily seen that the branching 
fraction can be approximated as 



B(Z^pie)~ (^) 2 ~10- 6 



(3.44) 



This is in agreement with its experimental value. Now we shall discuss these 
suppression in ji — > 3e decay. 
CASE-1: 

First discuss the vector-axial-vector form factor contribution at V^ite ver- 
tex to fi — * 3e, i.e. instead of the r y a we have 7° (fy + Ja 75)- Therefore the 
amplitude becomes 

, 9Vu e 9v ee 



A(/jl^ 3e) = u e (k 3 )>y a (f v + f A lb) u il (p)u e (k 2 )-f a v e (ki) 
= u e (k^ a (f v + f A >y 5 ) u^(p)u e (k 2 )'j a v e (ki) 



My-S 

9V„e9Ve e 
My 



Using the same procedure as for the previous cases \A\ , we get 



|A| 2 = 32 |/y| 2 + |/ A | 2 [(k 1 .p)(k 2 .k 3 ) + (k 1 .k 3 )(k 2 .p)] 



9V„e9Ve, 

, My 



(3.45) 



(3.46) 
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Using the spin averages, the above expression becomes 



E 14 

spins=±7£ 



\fv\ 2 + \f A \ 2 l(h.p) (fc 2 .fc 3 ) + (h.k 3 ) (k 2 . P )} 



, Ml 



(3.47) 



Solving the delta function using the technique described above, the decay 
width becomes 



T(fi^3e) = ^[\f v \ 2 + \f A \ 



The corresponding branching ratio becomes 



Ml 



7.2 ^2 



1927T 3 ' 



(3.48) 



B(^-3e)» \fv\ 2 + \f A \ 



9 v^9v ee i g w 

My ' Mi ' 



(3.49) 



CASE-2: 

Now consider the electric and magnetic form factor contribution to the 
ji — > 3e at Vi//e vertex. The amplitude can be written as 



A = M e (/c 3 ) 1 ^— (if M + fElb) (y av u^(p)u e (k 2 )-i a v e (k 1 ) ^ytt 1 



Mi 



(3.50) 



Using the same procedure as before the decay width for this reaction becomes 



r (n -> 3e) 



i i fgv^gve, 



2 80 v 



I/m| 2 + |/ e | 



M2, 1927T 3 



(3.51) 



The branching ratio for this case becomes 



B (fi -> 3e) 



80 



I/m| 2 + |/ 



m M ^ 9V»e9Vee / 9\V 

4 / iri 



(3.52) 



Hence, the branching ratios of the three body decays are suppressed by a 



factor ^k-. 
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Figure 3.1: /i — > 3e 



Figure 3.2: Z — > jie 
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Chapter 4 

Lepton Flavour Violation in 
LSND and Seesaw model of 
neutrino masses 

4.1 Introduction 

The question of whether or not neutrinos have a nonzero mass has remained 
one of the most tantalizing issue in the present day physics. In the Standard 
Model of electroweak theory neutrinos are considerd to be massless. But 
there is no compelling theoratical reason behind this assumption [1]. Hints 
that the neutrinos are massive particles comes from the observation of three 
anomalous efects, 

1. the LSND (Liquid Scentilator Neutrino Detector) excess [2, 3, 4], 

2. the atmospheric anomaly [5, 6, 7] and 
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3. the solar neutrino deficit [8, 9, 10, 11]. 

In particular, the atmospheric results are the most convincing ones. All 
the three effects can be naturally explained in terms of neutrino flavour os- 
cillations, which will occur when neutrino propagate through space, if there 
masses are non-degnerate and the weak and mass eigenstates are mixed. 

4.1.1 The neutrino flavour oscillation 

However, in order to explain all three experimentally observed effects in terms 
of neutrino flavour oscillation, one is forced to invoke additional sterile neu- 
trino states [12] to accommodate the very different frequencies of oscillations, 
given by three different mass squared Am 2 's, indicated by three different ef- 
fects. The existance of such neutrino is currently unresolved problem and 
clearly demonstrates that the neutrino sector is not fully understood. So 
due to the unappealing theoretical feature of a light sterile neutrino, it is 
interesting to look for alternatives that could explain the LSND excess with 
the known three light neutrinos. 

From a phenomenological point of view, we recall that the neutrino flavour 
oscillation hypothesis predicts a well defined dependence of phenomena as a 
function of neutrino energy, characterized by the so called L/E behaviour, 
where L is distance between source and detector and E is neutrino energy. 
So far, no experiment has conclusively demonstrated such a L/E dependence 
of the anomalous effect, with may be expectation of SuperKamiokande data 
which favours the dependence a LE n where n = — 1 [13]. 
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In such an unclear sitution, it is possible to envisage "non-flavour oscil- 
lation" mechanism to explain part of neutrino data. 

Aside from theoretical argument against sterile neutrino we argue that, 
from a phenomenological point of view, the LSND effect is peculiar: it has a 
small probability, measured to be (2.5 ± 0.6 ± 0.4) x 10 -3 [13], in contrast to 
the solar and atmospheric neutrino anomalies, which are large. Hence, LSND 
is a natural candidate for an interpretation involving a different physics than 
in atmospheric and solar neutrino flavour oscillation. 

4.2 LSND puzzle 

We recall that the LSND effect was first reported as an excess of z/ e 's in 
the flux from the /z + decay at rest (DAR) process. The neutrino beam is 
obtained with 800 kinetic energy protons hitting a series of targets, producing 
secondary pions. Most of 7r + come to rest and decay through the sequence 
7r + — > /U + z/ M ; followed by fi + — > e + u e u ai supplying the experiment with the 
z/^'s with a maximum energy of 52.8 MeV. The intrisnic contamination of 
z/ e 's coming from the symmetrical decay chain starting with n~ is estimated 
to be small since most negatively charged mesons are captured before they 
decay. 

The excess of z/ e 's, explained in terms of neutrino flavour transitions of 
the type u u — > u e , occur via the reactions: 

fi + — >■ e + u e u u 

- vacuum _ 
Vu ► V e 
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v e p — > e + n 

i.e. anti-neutrinos are produced by the /z + — > e+z/gZ/^ and detected by u e p — > 
e + n. 

There is another evidence in favour of the neutrino flavour oscillation 
which was reported in the Decay In Flight (DIF). But due to its lower statis- 
tical significance, we concentrate on the hint from stopped muon, and ignore 
the DIF result. 

The LSND claim is contradicted by the latest KARMEN2 results [14], 
however the experimental sensitivity is marginal to conclusively exclude or 
confirm completely the LSND excess. A new experiment, MiniBOONE [15] 
would confornt the flavour oscillation hypothesis with a very high statistical 
accuracy. A negative result from MiniBOONE experiment would indicate 
that the neutrino flavour oscillation is not the correct hypothesis to explain 
the excess seen in LSND. It would however not contradict other possible 
non-flavour-oscillation interpretation of the effect. In particular, LFV decays 
would play a role in the interpretation of LSND excess. A neutrino factory 
is an ideal mechaine to probe such anomalous decays of muon. 

Th aim of this work is to investigate another approch. We assume that 
the three neutrinos interact through Lepton Flavour Violation interactions, 
which are forbidden in the Standard Model. This is an attractive possibility, 
because various extensions of the Standard Model which predict the neu- 
trino masses also give rise to such new interactions. These interactions can 
affect the LSND excess. The motivation for the LFV is that the branching 
probability for LSND is too small ((2.5 ± 0.6 ± 0.4) x 10" 3 )) compared to 
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atmospheric and solar anomalies. We analyze the consequence of small Lep- 
ton Flavour Violating interactions to explain the LSND excess, and check 
whether this scenario will be feasible. 

4.2.1 Basic consideration and calculation 

Let us assume that a vector boson Vi (here Vi could be either a fundamental 
state, such as the Z° or the quark- antiquark bound state such as the 0, J/ip, 
or T) couples to /i ± e ± . If it couples also to e + e~, as all the states above 
do, then by unitarity its exchange contributes also to /i — > 3e. Let us define 
the effctive couplins between vector boson Vi and / u ± e =F as gv^, and the 
corresponding effective Lagrangian can be written as 

£e ff = ~gv^ a eV a + H.c. (4.1) 

Here, the reaction under consideration is 

therefore amoung these vector bosons only the Z° will contribute. / is any 
of the three known leptons, i.e. e, /i, or r. The coupling at the LFV vertex 
remains the same, but at the neutrino vertex it will be changed and is ob- 
tained by the Standard Model. We also assume that the LFV decay proceds 
through a similar diagram as the Standard Model muon decay, however with 
interchanged neutrino flavour as shown in fig. 4.3.1. 

Using Feynman rules, the corresponding amplitude can be written as 

v{pi)i a v{p2)u{pA)^ (1 - 75) V(p 3 ) 
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A (fi — ► evivi) 



Ml — c 



= f|r!f lHPlh a v(P2)u(P4ha(l-l 5 )v( P3 )} (4.2) 

Here gy u p is the effective coupling of the Z- boson to vv. Also for low energy, 
we have to neglect s in comparison to the M z . Therefore Eq. (4.2) takes the 
form 

A (/i -> evin) = 9Z ^ VV [v{pi)l a v{p 2 )u{p A )'-i a (1 - 75) v(p z )\ (4.3) 

z 

2 

Now as usual we have to calculate \A\ , which is 

2 



,2 J_ / 9Z^9Zvv 

16 V MI / 
xTr {M« (1 - 7s) M/3 (1 - 7s)} • (4.4) 

Solving these traces, the above Eq. becomes 

E l^l 2 = (n^) 4 [(PiP*) (PaPs) + (pips) (P2P4)] . (4.5) 

Then the partial decay rate for the muon can be written as 
1 1 d 3 k 3 fd^d 3 ^ 



dT 



mE J-^ S |j4| 2 (2it)V(pi +pi-q), (4.6) 
where 

<? = P3 + Pa. 

Using the same technique as we have used in the previous chapter, we got 
the total decay width as follows 

/ ~ \ 2 s 

w - \ / 9z^9zw \ m v „s 

V ^^ evm) = \-i^r) i92^- (4 - 7) 
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Comparing the above contribution to the \i — > eviv\ process to that of 
the ordionary muon decay, /x — > euP, which proceeds via W exchange and 
(almost) identical kinematics, gives the relation 



r (jj, -> eviv x 



Z—exch. 



r(/i evv) 
B(fi-> em) z _ exch . 



Z—exch. 



r (/i — > 3e) 



V — exch. 



r [/x -> 3e] v _ exch- 
Using the value of the r [// — > 3e] z _ exch , Eq. (4.8) becomes 



T(/i — > ez/z/) 
xB(^3e) M . (4.8) 



B (/x -> evivi) z _ exch . 
We know that 



0z* 



9Ze 



x B -> 3e) 



V— eircZi. 



(4.9) 



[^(/^^3e)] y _ e:cc ,.<10- 



-12 



Since 



and 



r (Z - ^) ~ <?f „ P M Z . 



So Eq. (4.8) becomes 



< 2 x 10 



-12 



(4.10) 



As it is seen from the above equation the branching ratio for the this lepton 
flavour violating decay is the same as that of the \x — > 3e. This is because 
of the fact that the reaction kinematic is same , and also in both the cases 
the external particles are assumed to be the massless. The factor of 2 is due 
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to the fact that the couplings between electron-positron is different then the 
two neutrinos. 

It is already mentioned that, in order to explain the LSND result, the 

effective New Physics coupling should satisfy 

r<v 2 



r = 



Gi 



= (2.5 ±0.6 ±0.4) x lO" 3 . 



(4.11) 



We define G V N to be the effective coupling of the anomalous muon decay. 
Thus, at the 90% C.L. we need [4] 



r > 1.6 x 10" 3 . 



(4.12) 



Now let's calculate this for the reaction discussed here. As we know that 

Gf 9w 



y/2 Mfa" 



(4.13) 



Thus following the Eq. (4.12), we can write the effective coupling as follows 

(4.14) 



G e J f gz^gzuu 



V2 



Now the ratio of the two couplings become 

G e / f gz^Qzuu 1 



So r becomes 



G F Ml G F 



G e J f |2 



(4.15) 



G F 

r (ji -> Ze) Z -exch. 1 9Zuu 



(4.16) 



As and using the value of the g\ vV and g 2 Ze+e - from the Standard Model, we 
get 



r ~ 10 



-12 



(4.17) 
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Thus, comparing with Eq. (4.12) which requires r > 1.6 x 10~ 3 , we learn that 
the LFV anomalous muon decays /i — > evxv x cannot significantly contribute 
to the LSND DAR result, at least when we use <7z +e _ as constrained by 
ji — > 3e in a model independent way. 

Thus the excess found by the LSND is not due to these Lepton Flavour 
Violating decays. So we cannot say that whether this excess is due to neutrino 
flavour oscillations or due to LFV decays which might occur through some 
exotic mechanism leading to /i + — > e + u e ui. Hence, we can say that it is still 
unresolved puzzle of the neutrino physics. 

4.3 Seesaw model of neutrino masses 

In the minimal standard gauge model of quarks and leptons, each of three 
known neutrinos (v e , u T ) appears only as a member of left-handed SU{2) 
lepton doublet 

^i = M) L , (4.18) 
and the Higgs sector contains only one scalar doublet 

$=(0 + ,0°). (4.19) 

As a result, neutrinos are massless in this model. Experimentally there is 
now a host of evidence for neutrino oscillations, and that is most naturally 
explained if neutrinos are massive and mix with each other. Theoretically 
there is no compelling reason for massless neutrinos, and any extension be- 
yond the minimal Standard Model often allows them to massive. There exists 
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already a vast literature on specific models of neutrino masses and mixing 
[16]. 

Here we will make the following simple observations. In the minimal Stan- 
dard Model of particle interactions, the massless neutrinos acquire naturally 
small Majorana masses through the effective dimension-five operator 

j (W° - - , (4.20) 

where A is a large effective mass scale, and $ = (0 + , 0°) is the usual Higgs 
doublet with the non zero expectation value, (0°) = v. All models of neu- 
trino mass and mixing (which have the same light particle content as at the 
minimal Standard Model) can be explained by the operator 

A~VW.r (4.21) 

Different models are merely different reasons for this operator. The interme- 
diate heavy particle in case of the operator defined in Eq. (4.20) is clearly a 
fermion singlet. Let's call it JVj and let its mass be mjy and its coupling to 
Vi be /j. Also we can identify j(v or simply fv as a Dirac mass m D linking 
V; L to Ni and the neutrino mass matrix is introduced so that 

m v = ^, (4.22) 

m N 

so that A = m^/f 2 in Eq. (4.20). This is, of course, just the well known 
canonical seesaw mechanism, with iVj identified as the right-handed neutrino 
with a large Majorana mass. Given that m v is at most of order leV and 
/ should be too small, the usual thinking is that m N has to be very large, 
i.e., m N ^> v. As such, this famous mechanism must be accepted on faith, 
because there cannot be any direct experimental test of its validity. 
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Consider now the possibility that there is no new physics beyond the 
TeV energy scale. This is an intriguing idea proposed recently in theories of 
large extra dimensions [17]. Instead of using an ingradient supplied by the 
large extra dimensions, a model has recently been proposed [19] to show how 
Eq. (4.22) may be realized naturally with of order 1 TeV in a simple 
extension of the Standard Model. This means that mjj should be small i.e., 
m D <C 10 2 GeV. If it comes from zeorth component of $ i.e., 0° as in the 
Standard Model, that would not be natural; but as shown below, it will come 
instead from another with naturally small vacuum expectation value. This 
new realization of the seesaw mechanism will allow direct experimental test 
of its validity, as discussed below. 

Consider the minimal Standard Model with three lepton families: 



\ */ L 

where their transformations under the standard SU(3) C x SU(2) L x U(1) Y 
gauge group are indicated above. Now add three neutral fermion singlets 



but instead of assigning them the lepton number L — 1, so that they can 
pair up with the lepton doublet through the interaction Nr(vl(P ~ 
L = is assigned to forbid this Yukawa term. To complete the model, a new 
scalar doublet, called leptoquark [18] 




l m ~ (1,1,-1), 



(4.23) 



N iR ~ (1,1,0), 



(4.24) 




(4.25) 



is introduced with lepton number L — — 1. Hence the terms 




(4.26) 
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appear in the Lagrangian. The effective operator of Eq. (4.20) for neutrino 
mass is then replaced by one with 77 instead of 0, and, if (77°) = u is naturally 
small, the corresponding scale A will not have to be so large and Mj of Eq. 
(4.26) may indeed be of order 1 TeV. 

The Higgs potential of this model is given by 



V = m 2 &<& + mlrfr) + -Ai ($ f $) 2 + -A 2 (77*77 



(4.27) 



+A 3 (77*77) + A 4 (V??) (77*$) + A*L ($ f ?7 + J7 f $) , 



where the n\ 2 term breaks L explicitly but softly. Note that, given the 
particle content of this model, the /i 2 2 term is the only possible soft term 
which also breaks L. 

The equations of constraints for (0°) = v and (77°) = u are 



m\ + Ait; 2 + (A 3 + A 4 ) u 2 + n\ 2 u 







u 



m\ + \ 2 u 2 + (A 3 + A 4 ) v 2 + [i\ 2 v = 0. 



Consider the case 



ml < 0, mj> 0, 



1A2 



'2: 



then 



7J 2 ~ 



(4.28) 
(4.29) 



(4.30) 



(4.31) 



Ai' m 2 + (A 3 + A 4 )t; 2 ' 

Hence it is very clear that u may be very small compared to v(— 174 GeV). 

For example, if m 2 ~ ITeV, \fi 2 12 \ ~ 10Ge\/ 2 , then u ~ lMeV. The rela- 
tive smallness of |/i 2 2 | may be attributed to the fact that it corrosponds to 
the explecit breaking of the lepton number in V of Eq. (4.27). The usual 
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argument here is that, if |//f 2 | were zero, then as seen from Eq. (4.31) the 
model's symmetry is increased, i.e., the lepton number would not be broken. 
Hence the assumption that it is small compared to \m\\ or m 2 is "natural." 
One thing is very clear from here that if \fJ,\ 2 \ were much smallar, then neu- 
trino masses would be too small to account for the present observation of the 
neutrino oscillations. 

The 6x6 mass matrix spaning [u e , u T , N±, N 2 , N 3 ] is now given by 



The mixing between v and iV is thus of the order fu/M, which will allow 
the physical N to decay through its small component of v to l ± W T . The 
effective mass matrix spaning the light neutrino is then 











felU fe2U / e3 M 



M v = 



fau f^ 2 u f^u 

f Tl u f r2 u f r3 u 

felU fau f Tl U M x 

fe2U f^U f r2 U M 2 

fesu fan f r3 u M 3 



(4.32) 




(4.33) 




ml + \ 3 v 2 + (A 2 - A 4 ) u 2 - n\ 2 u/v, 
m l + (^3 + A4) v 2 + \ 2 u 2 - n\ 2 u/v, 



(4.34) 



(4.35) 




(4.36) 




(4.37) 
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The m^o behaves very much like the ordionary Higgs boson. The new scalar 
particles of this model, i.e. h ± , A, and h® (all with mass~ m 2 ), as well as 
N iR , are now accessible to direct experimental discovery in future accelators. 

In summary, a new seesaw model of neutrino mass is proposed, where a 
second scalar doublet (77+, 77 ) with lepton number L = — 1 is added to the 
minimal Standard Model togather with three neutral right-handed fermion 
singlets iVj with lepton number L = 0. Thus iVj is allowed to have a Ma- 
jorana mass as well as interaction fijNm (vjLV° ~ ljLV + )- Hence m v is 
proportional to {if) 2 /rrii and, if (77°) <C (0°) , m N may be of the order 1 
TeV and be observable experimently. This is accomplished with the Higgs 
potential of Eq. (4.27), where L is broken explicitly and uniquely with the 
soft term &rj + 

As the Lepton Flavor Violation is discussed in all models of neutrino 
mass. It is argued that in this model, there is no LFV at tree level for charged 
leptons. However, it does occur in one loop throug 77 and N exchange [19]. 

The aim of this work is to test this model for the lepton flavor violating 
decays, i.e. /i — > 3e and /i + — > e + uiUi in one loop process. We have to 
calculate the branching ratio of \i — > 3e and compare it to its experimenal 
bounds. 

4.3.1 Basic consideration and calculation 

First of all we will discuss the reaction /i — > 3e. The corresponding box 
digrams can be shown in fig. 4.3.1. The couplings at each vertex can be 
taken from Eq. (4.32). Instead of 7 materices we have unity at each vertex, 
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because each Higgs boson is a scalar. 

Using Feynman rules, the amplitude can be written as follows [20]: 



1 t (/i - 3e) = 2^(f; i f:j:je i ) J ^«(p0 



where 



xm(j) 3 ) 



Jk — rrii 



v(p 4 ) 



jk — mi 

2 



k 2 — m\ 



xu(p 3 ) [/; + nn] v(p 4 ) 
2 



A; 2 — m : 



1 2 



(2vr) / 



/c 2 — m 2 



+rrii [v(p 1 )-f a v(p 2 )u(p 3 )v(p 4 )I a + v(p 1 )v(p 2 )u(p 3 )-f l3 v(p4)I 13 ] 

+771- {^(piM^MPsM^K}] , 



ja/3 

7 



2^ 2 



(A; 2 — mfj) (k 2 — m 2 ) 
k a 

(k 2 — ml) (k 2 — m 2 ) 
1 



(4.38) 

(4.39) 
(4.40) 
(4.41) 



(k 2 — mjjf (k 2 — m 2 ) 2 
It is assumed that the loop momenta is very high, i.e. k — > oo. So there is 
no change in momenta at the vertices, and also the dominating integral is 
the first one. Before calculating the |T| 2 , we have to do the loop integration. 
To do this we use the technique known as Feynman parameterization [21]. 



1 yi 

a 2 b 2 Jo 



dz- 



*(!-*) 
[b + (a- b)z 



4 ' 



(4.42) 
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For our above mentioned case a and b are 



k — rri; 



m h . 



(4.43) 



So 



[b + (a — b)z\ = k 2 — ml + {m 2 h — m 2 ^j z 
= k 2 + m 2 h [(1 — Xi)z- 1] 
= k 2 + s, (4.44) 

where X{ = (j^j and s = m\ [(1 — Xi) z — 1]. Then, Eq. (4.39) becomes 



J a/3 = xz(l-z) J d 4 k 

= 6 / dzx^(l-z) x 
Jo 

= 6 / (Izxz(l-z)x 
Jo 



[fc 2 + s] 



in 2 T(l) x # 
2r (4) x s 



a/3 



2 x 6m 2 [(l- Xi )z- 1]J ' 



(4.45) 



because 



/ 



m 2 r (n - 3) gf 



a/3 



n-3 



ra > 3. 



[k 2 + s + ie] n 2T(n) s" 

Therefore, after solving the z integration Eq. (4.45) takes the form 

jafs = ™ 2 g aP x \ l-x 2 + Xjhi (x 2 ) 
2m 2 [ 2 (,• 1 ) :! 

Using the same technique for Eq. (4.40) becomes 

k a 



(4.46) 



r 



<fk- 



-I 



(k 2 — m 2 h ) 2 (k 2 — m 2 ) 2 



d 4 k- 



k a 

(k 2 - s y 



= o, 



(4.47) 
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because 



Also define 



/ 



d 4 k 



k a 

[k 2 + s + te] n 



= n>3. 



/>* p 
fiiJeiJeiJei Si- 

Hence finally Eq. (4.38) becomes 



(4.48) 



T (fj, -> 3e) 



7T 



(27r) 4 m 2 

r 1 — x 2 + Xj In (a;. 2 



{ V (pi ) 7 Q V (p 2 ) « (P3 ) 7a W (p 4 ) } 



X 



I 2(x t - l) 3 
Now the square of the amplitude becomes 

' .2 \ 2 



(4.49) 



E m- 



7T 



spins 



X2n) 4 ml 



1 -x\ + Xi In (xf) 
2 (Xi - if 



x 



1 - x 2 + xi In (a; 2 ) 
2 - l) 3 



x2 [(P1P4) (P2P3) + (piPs) (P2P4)] • 



Define 



1 - x 2 + Xjln (xf) 



2 (a* - 1)' : 



Then the above result takes the form 



E m 



spins 
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The corresponding decay width becomes 



(4.50) 



r [// -> 3e] 



7T 



-1 2 



L(2vr) 4 m 2 j 



E^;-u.o!.u. 



x 



192 X 7T 3 ' 



(4.51) 
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and the branching ratio becomes 

r \ji — > 3e] 



£[//-> 3e] 



r [/i — > ez/z/ 
1 

(27r) 4 2m|. 



]T^*A(^fe) 



r 2 -> ez/] x 2 5 ' 



(4.52) 



Now let's calculate the other reaction, i.e. fj, + — > e + z//z//. 

The corresponding box diagrams can be shown in fig. 4.3.1. 
Using Feynman rules, the amplitude can be written as follows 



^(/x+^e+z^) = 2£(/^a)/^(p 



xv(p 2 )u(p 3 ) 



jk — ITLi 



v{p A ) 



]k — m, ( 
—i 



k 2 — m^+j 



—i 



k 2 — m 2 h0 



= 2 £ (fUfufl.fi) J 7^TI«(Pi) [/* + "a] v(p 2 ) 



xu{P3) [A + m] v(p A ) 



k 2 — m 2 h+ \k 2 — m 2 ) k 2 — m 2 h{) 



(4.53) 



Using the same assumption here that the loop moment is very high, i.e. 
k — > oo. The above integral can be solvd by using Feynman parameterization 
[22] 

— — = 6 dx dy— 

a z bc Jo Jo [a 

Therefore 



(1-x) 



ja/3 



I 



d A k 



(b — a) x + (c — b) y] 
k a k? 



-ITT 



Am 



(k 2 — ni 2 ) 2 (k 2 — m 2 h+ ) (k 2 
2 A( y x li ,x 2i )g a ' 3 , 



(4.54) 
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where 



and 



Define 



•"11 9 5 

mi 



%2i = 



m h° 2 ,A 

mi 



A(x U ,X2i) = 

J (x u ) = 



J (%) - J (x 2i ) 

Xli — %2i 
1 ~ 2 



3/ ^ In x ^2 
1 -a:w) ' (1 -x u y 



+ 



fiiJliJeiJli Si- 

Using these results Eq. (4.53) takes the form 



—m 



2m? (2tt) 4 . 

X {^(Pl)7^(P2)M(P3)7a^(P4)} • 



(4.55) 



(4.56) 



(4.57) 



Now, after calculating spins \T\ , the decay width becomes 



7T" 



2m? (2tt) 4 



to: 



32to 2 tt 2 



192 X 7T 3 



to: 



192 x 7T 3 ' 

(4.58) 



Then the branching ratio becomes 

1 2 



B 



e v x v x 



32mh 2 



Yl {&3 A ( X U, X 2 i)A(x lj ,X 2j )) 



2 x r 2 [W -> ei/] 
(4.59) 

In order to calculate the numerical value of the branching ratio for /i — > 3e, 
we will fix the value of TOj — 1 TeV and also assume that all the couplings 
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are of the order unity. As m\ ~ m|, therefore if m 2 ^> 100 GeV its effect on 
the radiative parameters 

1 X4V 2 



AS 



AT 



247T ml 

1 _1 Xjv 4 

967rs 2 c 2 Mf ml ' 



is negligible small and will not change the excellent experimental fit of the 
minimal Standard Model. 

To have a feeling about the branching ratios let us take mn = 900 GeV, 
the branching ratio for /i — > 3e becomes 

B = 3.6 x 10" 6 

which is much higher then the experimental bound on the branching ratio of 
/i — * 3e, i.e. B (/i — > 3e) y _ ea;cfe < 10 -12 . For the mass smallar than 900 GeV 
the branching ratio grows higher and higher. Thus from here it is obvious 
that mh > 1 TeV. By varing the Higgs mass the corresponding branching 
ratio becomes 

m h = lOTeV B = 1.0 x 10~ 8 

m h = 30TeV B = 3.1 x 10" 11 . 
It can be eaisly seen that with the increase of Higgs mass, the branching ratio 
decreases. As the Higgs mass reaches to 50 TeV, branching ratio becomes 



B = 3.6 x 10 



-12 



which is comparable to the experimental value. Thus we have concluded 
that the Higgs mass is greater than 50 TeV to explain LFV decays in Seesaw 
model of neutrino masses [19]. 
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Conclusion 



The aim of this work has been to understand the physics beyond the Stan- 
dard Model by considering Lepton Flavour Violating (LFV) decay processes. 
Using the experimental bounds on the three body LFV decays, we have cal- 
culated the bounds on the corresponding two body decays. The Dynamical 
suppression of three boby LFV decays due to momentum dependent cou- 
plings have also been discussed in detail. 

The experimental bounds on the three body LFV decay, /j, — > 3e con- 
strained the coupling gz^ e , which is helpful to calculate the anamolous muon 
decays. The anamolous muon decay fi — > ev\i>i could not significantly con- 
tribute to the LSND DAR result. This has been concluded by comparing the 
effective coupling of anamolous muon decay with r > 1.6 x 10 -3 [4]. LFV 
decays at a loop order have also been discussed in Seesaw model of neutrino 
masses which involve right handed singlet neutrinos of mass of order TeV 
[19]. The main purpose of the model is that there is no new physics beyond 
1 TeV. It has been shown that even if one keeps m N at 1 TeV the present 
bounds on /i — > 3e requires that new Higgs bosons, necessary in this model 
should have their mass larger than 50 TeV. Thus it is necessary to have a 
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new physics beyond the TeV energy scale defeating the origional motivation 
of the model. 
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Figure 4.1: Anamolous muon decay 



Figure 4.2: Box diagram for // — > 3e 



Figure 4.3: Box diagram for anamolous muon decay 
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